    % section 6.4.1 'robustness' / GKV 
        % this file shows the differences in CE for GHHW and GLTHI 
        % using Epstein Zin preferences

    clear;
    load('./results/allstates_main.mat')
    load('./results/model_parameters_main.mat')
    load('./results/capHtmod_c.mat');
    iX = simulation_output.iX;
    addpath functions
    replicate = 1;

    % SHI (aka GKV) starting distribution 
    shiprobs = importdata('../../data/processed/HealthProcess/acgdist25shi.csv');  
    val=shiprobs.data;
    % add mortality --
    capHtmod_shi = [val 0];
    % eliminate 7 and distribute equally
    capHtmod_shi(1,7)=0;
    capHtmod_shi(1,:)=capHtmod_shi(1,:)/sum(capHtmod_shi(1,:));
    
    % 90% SHI 10% PHI 
    capHtmod_mix = 0.9*capHtmod_shi+0.1*capHtmod_c;

    % SHI/GKV distribution

    if replicate ==1
        model_parameters.outputfile ='./results/allstates_shi.mat';
        model_parameters.capHtmod= capHtmod_shi;
        load('./results/iX_submission/iX_shi.mat');
        model_output_shi = solve_main(model_parameters,iX);
        writetable(model_output_shi.maintable,'../../figsandtabs/table_gkv.xlsx')

        model_parameters.outputfile ='./results/allstates_mix.mat';
        model_parameters.capHtmod= capHtmod_mix;
        load('./results/iX_submission/iX_mix.mat');
        model_output_mix = solve_main(model_parameters,iX);
        writetable(model_output_mix.maintable,'../../figsandtabs/table_mix.xlsx')
    else
        model_parameters.outputfile ='./results/allstates_shi.mat';
        model_parameters.capHtmod= capHtmod_shi;
        model_output_shi = solve_main(model_parameters);
        iX = model_output_shi.iX;
        save('./results/iX_out/iX_shi.mat','iX');
        writetable(model_output_shi.maintable,'../../figsandtabs/table_gkv.xlsx')

        model_parameters.outputfile ='./results/allstates_mix.mat';
        model_parameters.capHtmod= capHtmod_mix;
        model_output_mix = solve_main(model_parameters);
        iX = model_output_mix.iX;
        save('./results/iX_out/iX_mix.mat','iX');
        writetable(model_output_mix.maintable,'../../figsandtabs/table_mix.xlsx')
    end

    

